Dynamical Mean-Field Theory within an Augmented Plane- Wave Framework: 
Assessing Electronic Correlations in the Iron Pnictide LaFeAsO 
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We present an approach that combines the local density approximation (LDA) and the dynamical 
mean-field theory (DMFT) in the framework of the full-potential linear augmented plane waves 
(FLAPW) method. Wannier-like functions for the correlated shell are constructed by projecting 
local orbitals onto a set of Bloch eigenstates located within a certain energy window. The screened 
Coulomb interaction and Hund's coupling are calculated from a first-principle constrained RPA 
scheme. We apply this LDA-I-DMFT implementation, in conjunction with a continuous-time quan- 
tum Monte-Carlo algorithm, to the study of electronic correlations in LaFeAsO. Our findings sup- 
port the physical picture of a metal with intermediate correlations. The average value of the mass 
renormalization of the Fe 3d bands is about 1.6, in reasonable agreement with the picture inferred 
from photoemission experiments. The discrepancies between different LDA+DMFT calculations 
(all technically correct) which have been reported in the literature are shown to have two causes: 
i) the specific value of the interaction parameters used in these calculations and ii) the degree of 
localization of the Wannier orbitals chosen to represent the Fe 3d states, to which many-body terms 
are applied. The latter is a fundamental issue in the application of many-body calculations, such 
as DMFT, in a realistic setting. We provide strong evidence that the DMFT approximation is 
more accurate and more straightforward to implement when well-localized orbitals are constructed 
from a large energy window encompassing Fe-3d, As-4p and 0-2p, and point out several difficulties 
associated with the use of extended Wannier functions associated with the low-energy iron bands. 
Some of these issues have important physical consequences, regarding in particular the sensitivity 
to the Hund's coupling. 



PACS numbers; 71.15.Mb, 71.10.Fd, 71.20.Be, 74.70.-b 



I. INTRODUCTION 

This article has two purposes. The first one is to 
present a new implementation of dynamical mean-field 
theory (DMFT) within electronic structure calculation 
methods. This implementation is based on a highly pre- 
cise full-potential linear augmented plane wave method 
(FLAPW), as implemented in the Wien2k electronic 
structure codoi. The second purpose of this article is 
to report on DMFT calculations for the iron oxypnic- 
tide LaFeAsO, the parent compound of the Till '-family 
of recently discovered iron-based superconductors. The 
strength of electronic correlations in these materials is an 
important issue, which has been a subject of debate in 
the literature^i^iii^i^. 

The combination of dynamical mean-field theory with 
density-functional theory in the local density approxima- 
tion (LDA-I-DMFT) provides a powerful framework for 
the quantitative description of electronic correlations in 
a realistic setting. A number of materials have been in- 
vestigated in this framework over the past decade, such 
as transition metals and transition-metal oxides, rare- 
earth and actinide compounds, and organic conductors. 



These examples testify to the progress in our understand- 
ing of the key physical phenomena associated with the 
competition between the localized and itinerant charac- 
ters of electrons belonging to different orbitals (see e.g. 
Refs. SiiOinilll for reviews). 

In the past few years, a new generation of 
LDA+DMFT implementations have been put for- 
war d-'^'^'-'^'^d^i^^d'^'i-'^^d^ , which emphasize the use of Wan- 
nier functions as a natural bridge between the band- 
structure and the real-space description of the solid in 
terms of orbitals. These functions span the subset of or- 
bitals which are treated within the many-body DMFT 
framework. In this article, we present an implemen- 
tation of LDA-I-DMFT within the full potential lin- 
ear augmented plane wave (FLAPW) framework, us- 
ing atomic orbitals that are promoted to Wannier func- 
tions by a truncated expansion over Bloch functions fol- 
lowed by an orthonormalization procedure. This is a 
simpler alternative to the previous implementation of 
DMFT within FLAPWii, which constructed the Wan- 
nier functions following the prescription of maximal local- 
isation22i2i. The choice of FLAPW is motivated by the 
high level of accuracy of this all-electron, full-potential 
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method. In the present work, we use the Wien2k elec- 
tronic structure package^, and we have constructed an 
interface to it that allows for the construction of Wannier- 
like functions used in DMFT. Our implementation is 
described in detail in Sec. [TTl As a benchmark, we 
perform calculations on a test material. SrVOs, which 
are presented in Appendix |A] and compared to previ- 
ously published results for this materiali^iii'2^'^'2ii2^i2£. 
Throughout this article, many-body effects are treated 
in the DMFT framework using the recently developed 
continuous-time strong-coupling Quantum Monte Carlo 
algorithm of P.Werner and coworker o^^'^^ . Because very 
low temperatures can be reached, and very high accuracy 
can be obtained at low-frequency, this algorithm repre- 
sents a major computational advance in the field. 

In Sec. mil we address the issue of electronic cor- 
relations in LaFeAsO. The DMFT calculations which 
have been published soon after the experimental discov- 
ery of superconductivity in the iron oxypnictides have 
provided seemingly contradictory answers to this ques- 
tion. In Refs. [201 J K. Haule and G. Kotliar proposed 
that LaFeAsO is a strongly correlated metal, rather close 
to the Mott metal-insulator transition, and character- 
ized by a reduced value of the quasiparticle coherence 
scale, resulting in bad metallic behavior. In contrast, in 
Refs. [3001, V.Anisimov and coworkers proposed that 
these materials are in a weak to intermediate regime of 
correlations. 

Our LDA+DMFT calculations for LaFeAsO support 
the physical picture of a metal with intermediate corre- 
lations. The average value of the mass renormalization 
of the Fe 3c? bands is about 1.6, in reasonable agree- 
ment with the picture inferred from photoemission ex- 
periments. We also find that there is no technical in- 
consistency between different DMFT results reported for 
LaFeAsO before. We show that the discrepancies in the 
literature are due to two causes: i) the specific value of 
the interaction parameters used in these calculations and 
ii) the degree of localization of the Wannier orbitals cho- 
sen to represent the Fe 3d states, to which many-body 
terms arc applied. 

In Sec. IIIH we perform detailed comparisons between 
LDA-fDMFT calculations performed with different de- 
gree of localization of the correlated orbitals, associated 
with different choices of energy windows for the Wannier 
construction (and accordingly, different degrees of screen- 
ing of the interaction parameters). We point out sev- 
eral difhculties associated with the use of more extended 
Wannier functions associated with the low-energy iron 
bands only. Some of these issues have important physi- 
cal consequences, in particular regarding the sensitivity 
to the Hund's coupling. 

This article ends with several appendices, reporting 
on more detailed aspects or technical issues. Appendix 
A is devoted to a benchmark of our implementation on 
a "classical" test compound, SrVOs. Appendix B de- 
tails some technical issues associated with the projection 
scheme used to display partial spectral functions with a 



given orbital character. Appendix C discusses the influ- 
ence of spin-flip and pair-hopping terms on the degree of 
correlations, on the basis of model calculations. We con- 
clude that while these terms are indeed important close 
to the Mott transition, they can safely be neglected in 
the regime of correlations relevant to LaFeAsO. 



II. THEORETICAL FRAMEWORK 

A. Implementation of LDA+DMFT in the APW 
framework 

1. LDA+DMFT in the basis of Block waves 

To make this article self-contained and in order to de- 
fine the main notations, this subsection begins by briefly 
reviewing some essential aspects of the LDA-I-DMFT 
framework. The presentation is close to that of Refs. (T3 . 
[igj . where additional details can be found. 

Dynamical mean-field theory is a quantitative method 
for handling electron correlations, which can be described 
as an "effective atom" approach. The self-energy in the 
solid is approximated by that of a local model, a general- 
ized Anderson impurity model describing a specific set of 
atomic-like orbitals coupled to a self-consistent environ- 
ment. The self-consistency requirement is that the local 
on-site Green's function of the solid, calculated using this 
local self-energy, must coincide with the Green's function 
of the effective impurity model. 

In order to formulate the local effective atom prob- 
lem, a set of (orthonormal) local orbitals IXm"^), and cor- 
responding Wannier-like functions Iw^^), must be con- 
structed. These Wannier functions span the "correlated" 
subspace C of the full Hilbert space, in which many-body 
correlations (beyond LDA) are taken into account. This 
set of orbitals spanning the correlated subspace must be 
clearly distinguished from the full basis set of the prob- 
lem, in which the Green's function of the solid can be 
expressed. Obviously, the basis set spans a much larger 
Hilbert space, involving all relevant electronic shells. 

Below, we discuss in details how the Iw^^*^) are con- 
structed from the local orbitals Ixm'^)- The index m is 
an orbital index within the correlated subset, a denotes 
the atom in the unit cell, and a is the spin degree of free- 
dom. Projections of quantities of interest on the subset 
C are done using the projection operator 

P"-(k)=5]K-)(zz;^-|. (1) 

The effective impurity model is then constructed for 
the correlated subset C. It is defined by the Green's 
function of the effective environment, Q^rn'i^^ri) and 
Hubbard-Kanamori interaction parameters Umm'm"m"' ■ 
By solving this model in a suitably chosen manner one 
obtains the impurity Green's function G'^™?{iuJn) as well 
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as the impurity self-energy 
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(2) 



For the formulation of the self-consistency condition 
relating the lattice Green's function of the solid to the 
impurity model, it is convenient to choose the Bloch basis 
iV'ki/) ^^'^ complete basis set of the problem, since it is 
a natural output of any electronic structure calculation. 
The (inverse) Green's function of the solid expressed in 
this basis set is given by: 

^^(k, iuJn)^l' = [ii^n + - ^l.u)^uy' - Ki^' ^^n), (3) 

where are the Kohn Sham eigenvalues and 

S^j^, (k, is the approximation to the self-energy ob- 
tained by the solution of the DMFT impurity problem. It 
is obtained by " upfolding" the impurity local self-energy 



a, mm' 

(4) 

where P^-^(k) = i/) ^I'e the matrix elements of 

the projection operator, Eq. U]) and 



(5) 



Here, S^™/^ is the impurity self-energy, Eq. ([2]), ex- 
pressed in the local orbitals. and is a double- 
counting correction, which will be discussed in Sect. Ill 5] 
The local Green's function is obtained by projecting 
the lattice Green's function to the set of correlated or- 
bitals m of the correlated atom a and summing over the 
full Brillouin zone, 

= E nr(k)GL.(k,zc.„)P;„7(k). (6) 



Note that the local quantities C!^^,{iujn) and 
AE^™^(ia;„) carry also an index a, which we sup- 
pressed for better readability. 

The self-consistency condition of DMFT imposes that 
the local Green's function, Eq. ^ must coincide with the 
one obtained from the effective impurity problem. 
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(7) 



This equation implies that the Green's function of the ef- 
fective environment. Go, must be self-consistently related 
to the self-energy of the impurity model through: 
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(8) 



where the dependence of Gioc on Sjmp is specified by 
Eqs. (|3l6p . In practice, the DMFT equations arc solved 
iterativcly: starting from an initial Qq, the impurity 
model is solved for Simp, and a new Gq is constructed 
from ([5]). The cycle is repeated until convergence is 
reached. 



In order to construct the set of Wannier functions, we 
start from a set of local atomic-like orbitals IXm"^) defined 
in the unit cell. These orbitals can be expanded over the 
full Bloch basis-set as: 



IXkm) 



(9) 



This expansion is then truncated by choosing an energy 
window W, and restricting the sum to those Bloch states 
with Kohn-Sham energies e^i, within W. The number of 
bands included in W will in general depend on k and a. 
We thus define the modified orbitals (which do not form 
an orthonormal set because of the truncation): 



(10) 



Let us denote the matrix elements of the projection op- 
erator for this subset as 



(11) 



The matrix P^^f (k) is not unitary, except when the sum 
in (|10p is carried over all Bloch bands. It is also important 
to note that the matrices P"''^ are in general non-square 
matrices. They reduce to square matrices only in the 
case when the number of Kohn-Sham bands contained in 
the chosen window equals at every k-point the number 
of correlated local orbitals to be constructed. 

The orbitals |Xkm) '^^^^ orthonormalized, giving a 
set of Wannier- like functions: 



km 



) = E^: 



m.m' 



(12) 



where , = {0(k,a)-V2}-"^, and 0-':„,(k,a) 
(Xkml Xknv) Ihe overlap matrix elements. 
The overlap 0"j" j/ (k, a) finally reads 



o:i'(k, ^) = E (k)c:;"(k), (13) 

w 

while the orthonormalized projectors are then written as 
P-,-(k) = E {[O(k,a)]-i/n"'\o;(k). (14) 



2. Augmented plane waves 

In this work, the Bloch basis iV'k^) are expanded in 
augmented plane waves (APW/LAPW), which will be 
briefly described in the following. As was first pointed 
out by Slater—, near atomic nuclei the crystalline poten- 
tial in a solid is similar to that of a single atom, while 
in the region between nuclei (in the interstitial) the po- 
tential is rather smooth and weakly-varying. Hence, one 
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may introduce a set of basis functions, augmented plane 
waves, adapted to this general shape of the potential. 
First the crystal space is divided into non-overlapping 
mufhn-tin (MT) spheres centered at the atomic sites 
and the interstitial region in between. In the intersti- 
tial region(/) the APW ct>^{r) is simply the correspond- 
ing plane wave for given reciprocal lattice vector G and 
crystal momentum k: 

0^(r) = 47e^(''+«> re I, (15) 

where V is the unit cell volume. This plane wave is 
augmented inside each of the MT-spheres by a combi- 
nation of the radial solutions of the Schrodinger equa- 
tion in such way that the resulting APW is continuous 
at the sphere boundary. The APW arc then employed 



to expand the Kohn-Sham (KS) eigenstates for 
the full Kohn-Sham (KS) potential, (without any shape 
approximation) . In the original formulation of the APW 
method the radial solutions expanding a KS eigenstate 
inside MT-sphcres had to be evaluated at the correspond- 
ing cigenenergy leading to an energy-dependent basis set 
and, hence, to a non-linear secular problem. In order to 
avoid this complication, linearized versions of the APW 
method have been proposed. There are two widely-used 
schemes for the APW linearization. In the first, the linear 
APW (LAPW) method^, the plane wave is augmented 
within MT-spheres by a combination of the radial solu- 
tions, evaluated at chosen linearization energies En, and 
their energy derivatives. The resulting linear augmented 
plane wave then reads: 



6G(r) 



1 „i(k+G)r 
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MT 



(16) 



where the index a = 1, ...,Nat runs over all Na atomic 
sites in the unit cell, the coefficients Aim and Bim are 
determined from the requirement for the linear APW to 
be continuous and diffcrentiable at the sphere boundary, 
r and f are the radial and angular parts of the position 
vector, respectively. The energy-independent basis set 
(|16p leads to a linear secular problem, however, compared 
to the energy-dependent APW, a larger number of the 
LAPW in the basis set is generally required to attain 
the same accuracy. In order to decrease the requirement 
for the number of APW another linearization scheme, 
APW-f^o2i has been proposed. The APW-l-Zo basis set 
consists of the augmented plane waves evaluated at a 
fixed energy Eu: 



1 ^»(k-|-G)r 

'^^(r) - <! EAi^^''ur{r,E-)Yiir) 

Ini 



r e I 

r G -RmT' 



(17) 

where the coefficient Aim are determined from the re- 
quirement for (f>^{r) to be continuous at the sphere 
boundary. To increase the variational freedom of the 
APW-|-^o basis set the fixed-energy APW (fT7|) are supple- 
mented for the physically important orbitals (with / < 3) 
by the local orbitals (lo) that are not matched to any 
plane wave in the interstitial and are defined only within 
the muffin tin spheres (r G R%it) 

</'t,a(r) = [Ar™°"r''(^> E^i) + Br^^urir, E^i)] Yiif) 

(18) 

with the coefficients Aim and Bim chosen from the re- 
quirement of zero value and slope for the local orbital at 



the sphere boundary. 

Additional local orbitals (usually abbreviated with the 
capital letters as LO, (pf'ma) '^^^ tie introduced to ac- 
count for semicore states. They have a similar form as 
Eq. with the redefined A?'^'^ and a second term with 



a coefficient CJ 



Q,LO 
Ira 



and the radial function evaluated at 



a corresponding energy E21 for the semicore band. The 
coefficient Bim is set to in the APW-|-^o framework. 

Generally, in the full potential augmented plane wave 
method the LAPW, APW+Zo and LO types of orbitals 
can be employed simultaneously. The Kohn-Sham eigen- 
state is expanded in this mixed basis as: 



N, 

E 

i=l 



(19) 



where Nf, is the number of the orbitals in the basis set. 
The LDA-I-DMFT framework introduced in the present 
work can also be used in conjunction with any mixed 
APW+/0/LAPW/LO basis set. 



3. Local orbitals and Wannier functions in the APW basis 

Having defined the basis set we may now write down 
the expression for the Bloch eigenstate expanded in the 
general APW basis For r in the interstitial region, 

it reads: 



Npv 



G 



(20) 
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while for the region within the MT-spheres r G R%it ('^ ~ 
1, ...,Nat), we have: 



Npv 



k+G a,tT 



G 

Nlo 



(21) 



where A^'p^y is the total number of plane waves consid- 
ered in the interstitial which in turn is augmented inside 
each MT-sphere, Nio is the number of (j)]'^ ^ (r) orbitals of 
Eq. (jlSp and N^o the corresponding number of auxiliary 
orbitals for semicore states (pi^ ai'^)- 

In the framework of the APW method one has sev- 
eral choices for the 'initial' correlated orbitals \xZf)- 
Any suitable combination of the radial solution of the 
Schrodinger equation and its energy derivative for a given 
correlated shell {q:,^} can be employed, for example, the 
lo-orbital (HH]). In the present paper we simply chose 
the \xZfY^ as the solutions of the Schrodinger equation 
within the MT-sphere at the correspond- 

ing linearization energy En. 

Inserting \xZf) = ^"''^(^lO^m) ^'^'^ the expansion 



nj of the Bloch eigenstate in terms of APWs into 
Eqs. (jlOllip . and making use of the orthonormality of 
the radial solutions and their energy derivatives: 



\Eu)Yi, 



= Sll'mm' (22) 

- 0, (23) 



one obtains the following expression for the projection 
operator matrix clement: 

P^';(k) = {u'f %En)Yl, IV'L) 

^ (24) 

In this expression, the first term in the right hand side of 
p4)) is due to the contribution from the LAPW and/or 
APW-l-^o orbitals 



G 

Nia , Nlo 

+ 2^ '^lo ^Im + 2^ '^LO^lm 
"(0 = 1 nto = l 



(25) 



and the contribution due to the LO (semicore) orbitals 
that arises due to mutual non-orthogonality of the radial 
solutions of the Schrodinger equation for different ener- 
gies 

czM'^^'^) = 4'S^^r™^°o-:.™. (26) 



where Of^i,^, is the corresponding overlap: 

Or.„, = (^'^liJiO^l urf{Ei^,o)Yi)^Q (27) 



Then we orthonormalizc the obtained local orbitals to 
form a set of Wannier-like functions, Eq. The cor- 

responding projection operator matrix elements (j24p are 
orthonormalized accordingly using Eq. (|14p . 



B. Implementation and computational methods 

1. FL APW code 

For the electronic structure calculation we use the full 
potential APW+/o/LAPW code as implemented in the 
Wien2k package^. We have built an interface that con- 
structs the projectors to the correlated orbitals (P^'^(k)) 
out of the eigenstates produced by the Wien2k code, as 
described in Sect. Ill A 31 In order to obtain the local 
Green's function, the summation over momenta, Eq. ([6|), 
is done in the irreducible Brillouin zone (BZ) only, sup- 
plemented by a symmetrization procedure which is stan- 
dard in electronic structure calculations, 



BZ 



N, IBZ 



EA(k)=EE^^A(k)0^ 



(28) 



s=l k 



where A(k) is any k-dependent matrix, Ng the number 
of symmetry operations and Og the symmetrization ma- 
trices. Furthermore, we construct the local orbitals in 
the local coordinate system of the corresponding atom. 
This means that the equivalent atoms in the unit cell 
for which the DMFT should be applied, e.g. the two Fe 
atoms in the oxypnictidcs, are exactly the same and the 
impurity problem has to be solved only once. Afterward, 
the Green's function and self-energies are put back to 
the global coordinate system of the crystal in which the 
Bloch Green's function, Eq. is formulated. 
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2. Continuous-time quantum Monte-Carlo 



For the solution of the impurity problem we use the 
strong-coupling version of the continuous-time quantum 
Monte Carlo method (CTQMC)^'^. It is based on a 
hybridization expansion and has proved to be a very ef- 
ficient solver for quantum impurity models in the weak 
and strong correlation regime. It allows us to address 
room temperature (/3 = 1/kT « 40 eV^^) without prob- 
lems. In our calculations, we used typically around 
5 • 10^ Monte-Carlo sweeps and 1000 k-points in the ir- 
reducible BZ. Since the CTQMC solver computes the 
Green's function on the imaginary-time axis, an ana- 
lytic continuation is needed in order to obtain results 
on the real-frequency axis. Here, we choose to perform a 
continuation of the impurity self-energy using a stochas- 
tic version of the Maximum Entropy method^^ yield- 
ing real an imaginary parts of the retarded self-energy 
ReS(a; -f iO"*"), Iml](a; -t- iO+) which can be inserted into 
Eq. ^ in order to obtain the lattice spectral function 
and density of states. 



the matrices via 

[/ = — V C/"^* , 



N 



J = U- 



1 ^ 

V 



N{N- 1) 



(30) 
(31) 



As mentioned above, the LDA-I-DMFT scheme (as the 
LDA-I-U one) involves a double-counting correction SjJ^^/ 
in Eqs. (|4I5|) . Indeed, on-site Coulomb interactions are 
already treated on mean-field level in LDA. Several forms 
of the double-counting correction term have been pro- 
posed and investigatedi^^i^i^^ In this work we follow 
Ref. m and use the following double-counting correction: 



(32) 



where U is the average Coulomb interaction, J the 
Hund's rule coupling. iV^ the spin-resolved occupancy 
of the correlated orbitals. and Nc = A^J + N^. We com- 
pared the results obtained with Eq. ((32|) also with the 
double-counting correction given in Ref. [s^ which gave 
very similar results. 



3. Many-body interactions 

The CTQMC strong-coupling algorithm can deal with 
the full rotationally invariant form of the interaction 
hamiltonian^S. However, most calculations presented 
in this article will consider only the Ising terms of the 
Hund's coupling, yielding the interaction Hamiltonian: 



Hint — 2 / J U jj^r^iTlnia^m' a 

^ (29) 

mm' 

with U^^i and U^^, the reduced interaction matrices 
for equal and opposite spins, respectively. This enables 
us to take advantage of a maximal amount of conserved 
quantum numbers and, hence, perform the CTQMC cal- 
culation without any truncation of the local basis. The 
effects of spin-flip and 'pair-hopping' terms in the Hund's 
interaction arc discussed in Appendix [C] 

In our approach, the interaction matrices are expressed 
in terms of the Slater integrals , F^, and F^, where for 
d-electrons these parameters are related to the Coulomb 
and Hund's coupling via U = F^ , J ^ {F^ + F'^)/14, 
and F"^ / F^ ~ 0.625i^ Using standard techniques the 
four-index [/-matrix is calculated, and the reduced inter- 
action matrices are then given by f/^m' ^ Umm'mm' and 
Umm' = Umm'mm' " C^mm'm'm- With the abovc defini- 
tions, the Coulomb parameters U and J are related to 



C. Choice of energy window, localization of 
Wannier functions, and screening 

The Wannier functions defining the correlated sub- 
space of orbitals for which a DMFT treatment is per- 
formed, arc constructed by truncating the expansion of 
the initial atomic-like local orbitals to a restricted en- 
ergy window W, as described above. The choice of this 
energy window is an important issue, which deserves fur- 
ther discussion. Indeed, it will determine the shape and 
the degree of localization of the resulting Wannier-like 
functions. 

Let us consider first the case of a rather small energy 
window, containing only those bands that have domi- 
nantly an orbital character which qualifies them as "cor- 
related" (e.g., the Fe-3d orbitals in LaFeAsO or the V-Sd- 
t2g orbitals in SrVOa). In that case, the dimension of the 
Kohn-Sham Hamiltonian used in Eq. ([3]) coincides with 
that of the correlated subspace C (i.e. with the number 
of orbitals involved in the effective impurity model, for a 
single correlated atom per cell). The Wannier- like func- 
tions are then quite extended in real space, and resemble 
strongly the Wannier orbitals constructed within other 
schemes, such as the maximally localized Wannier con- 
struction of Ref. [l^, or the A^-th order muffin-tin basis 
set downfolded to that set of bands j^s^^i^ In such a situ- 
ation, hybridization of the correlated orbitals with states 
that lie outside the energy window is neglected at the 
DMFT level. Some information about the hybridization, 
e.g., of the d-states with ligand orbitals is of course taken 
into account through the leakage of the Wannier orbitals 
on neighboring hgand sites (see, e.g., Ref. |40[). 
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In contrast, if a larger energy window is chosen, it will 
in general contain states treated as correlated as well as 
states on which no Hubbard interactions are imposed. In 
this case, the dimension of the Kohn-Sham Hamiltonian 
used in the LDA+DMFT calculation of the local Green's 
function (j3l6|) exceeds the number of correlated orbitals 
involved in the effective impurity model. The Wannicr 
functions are more localized in space, and the information 
about the hybridization of the correlated orbitals with 
other states within this larger window is carried by the 
off-diagonal blocks of the Hamiltonian between correlated 
and uncorrclatcd states. 

An instructive case occurs when the correlated bands 
are well separated from the uncorrelated bands at all k- 
points, but the bands overlap in energy. This situation 
is realized for instance for the t2g bands in SrVOs that 
extend into the energy region of the eg bands. In order 
to strictly pick the three correlated t2g bands at each 
k-point, one would in that case have to introduce a k- 
dependent energy window. For a k-indepcndent energy 
window, one will in general have more than three bands 
at some k-points. corresponding to some Cg contribution 
in the chosen window. This can then be expected to 
result in slightly more localized orbitals. An example is 
given in Appendix \X[ 



D. Local Coulomb interactions, screening and 
constrained RPA calculations 

The choice of the energy window influences the value of 
the interaction parameters Umm'm"m"' in a crucial man- 
ner, which can be traced back to two main reasons. 

First; the interaction parameters are related to ma- 
trix elements of a screened interaction between the cho- 
sen Wannier functions. The more bands are included in 
their construction, the more localized they become and, 
hence, the matrix elements increase. Second, screening 
effects themselves affect the value of U . The more states 
are excluded from the screening process, the larger U be- 
comes. In what follows we distinguish carefully between 
these two effects. 

In the present work, we apply the present 
LDA-fDMFT implementation to one of the new 
high-Tc superconductors, LaFeAsO. For constructing 
the Wannier functions, we focus on an energy window 
that contains the 10 bands around the Fermi level with 
dominantly Fe-3(i character and also the bands coming 
from the p-bands of O and As, which are mainly located 
in the energy region [— 6,— 2]eV, resulting in a dpp- 
hamiltonian" . In addition, we performed calculations 
also for a smaller energy window containing only the 
Fc-3d bands, yielding a " d-hamiltonian" , as well as for a 
very large window including around 60 Bloch bands. 

The values of the Coulomb interactions U and J are 
calculated from the constrained random phase approxi- 
mation (cRPA) ^^i^^ , using the recently developed scheme 
for entangled band structures^. cRPA calculations for 



LaFeAsO have been performed before in Refs. 
Ref. \M the screened Coulomb parameters are obtained 
for three different situations: (i) by constructing Wan- 
nier functions from an energy window comprising the Fc- 
d bands only, and screening calculated excluding the Fe-c? 
channels only, (ii) by considering a larger window which 
also includes the As and 0-p states, so that the screening 
processes, for instance, from the As-p states to the Fe-c? 
ones are also excluded, (iii) and finally, a hybrid situa- 
tion (dubbed 'd-dpp' in Ref. IZsh . in which the Wannicr 
functions are calculated from an energy window includ- 
ing Fc-d, 0-p, and As-p, but only the Fe-d states are 
excluded from the screening. In other words, the screen- 
ing is calculated as in (i) and the Wannier functions as in 
(ii) . As discussed in Ref. \M this third option should be 
appropriate in a situation in which a full dpp-hamiltonian 
is used but Hubbard interactions are only applied to the 
d states. 

Here, we follow the same approach as in the "d- 
dpp" case of Ref. H^, but using the new disentanglement 
scheme of Ref. [isl. First, a partially screened Coulomb in- 
teraction Wr is constructed as follows: Wannier functions 
for the Fe-c? states are calculated, and a basis for the com- 
plementary subspace (containing in particular the lig- 
ands, but also higher lying f-statcs) is constructed. Based 
on the interpolating c?-band structure, the Pd polarization 
is computed, and the partially screened Coulomb inter- 
action Wr is obtained by screening the bare Coulomb 
interaction by all RPA screening processes except Pd- Fi- 
nally, according to the cRPA procedure the Hubbard U 
matrix is composed of the matrix elements of Wr in the 
basis of dpp- Wannier functions. As argued in Ref. \M 
this procedure is suitable for calculations that deal with 
the full c?pp-Hamiltonian in the many-body calculations 
while explicitly retaining Coulomb interactions on the d- 
submanifold only. In particular, we stress that to the ex- 
tent that our projection method produces Wannicr func- 
tions for the c?pp-window, the Hubbard U parameters are 
expressed in the same basis as the impurity quantities. 

It is important to note that we keep the screening 
channels in cRPA, i.e. P^, unchanged when the energy 
window in our calculation is varied. Thus, the differ- 
ent energy windows affect only the localization of the 
Wannicr functions, but not the screening process of the 
bare Coulomb interaction, and therefore the effective in- 
teractions are increasing with increasing energy window. 
Keeping the screening channels fixed is fully consistent 
with the fact that correlations are only included for the 
d electrons, but not for the ligand states. 

For our purposes, we calculate the average Coulomb 
interaction U and Hund parameter J from the matrices 
calculated by cRPA. With this U and J, the interaction 
matrices in the spherical symmetric approximation used 
in our calculation arc obtained as discussed above. As 
we will discuss in more detail below, the comparison of 
the resulting U^^^^^, and U^^^, with the cRPA matrices 
shows that for the dpp-liamiltonian, the approximation 
using atomic values for the ratios of Slater integrals F'^ 
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is well justified, whereas for the d-hamiltonian the cRPA 
matrices show strong orbital anisotropies. 



III. 



RESULTS FOR THE IRON OXYPNICTIDE 
LAFEASO 



Construction of the 



-Hamiltonian 



Let us start the discussion of correlation effects in 
LaFeAsO with our results for the dpp-hamiltonian, for 
which Wannier functions arc constructed from the en- 
ergy window W = [—5.5,2.5] eV. These Wannier func- 
tions are quite well localised. The corresponding Kohn- 
Sham hamiltonian contains 22 Bloch bands, correspond- 
ing to the 10 Fe-3(i bands, the 6 As-p and the 6 0-p 
bands. 

The local many-body interactions corresponding to 
this choice of Wannier functions are obtained from cRPA, 
as described in the previous section. They read: 



aa I 

mm' IcRPA 



/O.OO 1.61 1.55 2.26 2.26 \ 

1.61 0.00 2.50 1.82 1.82 ' 

1.55 2.50 0.00 1.70 1.70 

2.26 1.82 1.70 0.00 1.74 

V 2.26 1.82 1.70 1.74 0.00 / 



/ 3.77 
' 2.35 



IcRPA 



2.21 
2.71 
2.71 



2.35 
3.94 
2.87 
2.44 
2.44 



2.21 
2.87 
3.31 
2.29 
2.29 



2.71 
2.44 
2.29 
3.48 
2.29 



2.71 \ 
2.44 ' 
2.29 
2.29 
3.48 / 



The ordering of orbitals in those matrices is ^^2, d^2_y2, 



dxz-i d: 



yz- 



According to the conventions of the formulae Eqs. ((30)) . 
(Plj) . these matrices correspond to the values: U = 
2.69 eV and J = 0.79 cV. Using these values of U and J, 
we construct the spherically symmetric interaction ma- 
trices, 



U, 



(7(7 

mm' 



I 0.00 
1.49 
1.49 
2.30 

V 2.30 

/ 3.59 
2.19 
2.19 
2.73 

V 2.73 



1.49 
0.00 
2.57 
1.76 
1.76 

2.19 
3.59 
2.91 
2.37 
2.37 



1.49 
2.57 
0.00 
1.76 
1.76 

2.19 
2.91 
3.59 
2.37 
2.37 



2.30 
1.76 
1.76 
0.00 
1.76 



2.30 \ 
1.76 ' 
1.76 
1.76 
0.00 / 

2.73 \ 



2.73 

2.37 2.37 
2.37 2.37 
3.59 2.37 
2.37 3.59 / 



It is obvious, that the approximation of the cRPA ma- 
trices by using spherical symmetrisation is well justi- 
fied in this case, with the largest absolute deviation be- 
ing AJ7 K, 0.35 cV, corresponding to a relative error of 
around 0.09. The reason for this good agreement is that 
in the present case the Wannier functions arc already 
very close to atomic-like orbitals. We also checked that 
cRPA yields a significantly larger value of V for iron- 
oxide (FeO), as expected physically. 




FIG. 1: (Color online) Total DOS for LaOFeAs, dp-p Hamil- 
tonian. Black line: LDA DOS. Red line: LDA-fDMFT DOS. 




FIG. 2: (Color online) Real (top) and imaginary (bottom) 
part of the orbital dependent self-energy in the d-pp Hamilto- 
nian for U = 2.69, J = 0.79. 



B. LDA+DMFT Results {dpp hamiltonian) 

We carried out LDA+DMFT calculations for the dpp- 
hamiltonian using the above matrices at an inverse tem- 
perature = 40eV~^ (room temperature T = 300K), 
using the experimental crystal structure of LaFeAsO. In 
Fig. [T] we display the resulting total densities of states 
(DOS) together with the corresponding LDA DOS. The 
total densities of states were computed from the lattice 
Green's function, Eq. ([3]), traced over all ;^ e W and in- 
tegrated over BZ. In order to obtain the corresponding 
LDA densities of states EJ^^,, (k, i(jj„) in Eq. ^ was set to 
zero. 

One sees in Fig. [Uthat the LDA+DMFT DOS near the 
Fermi level displays characteristic features of a metal in 
an intermediate range of correlations. Both occupied and 
empty states are shifted towards the Fermi level due to 
the Fermi-liquid renormalizations. No high-energy fea- 
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FIG. 3: (Color online) DOS of Fe orbitals in LaOFeAs, dpp 
Hamiltonian. Color coding as in Fig. [T] 
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tures that would correspond to lower or upper Hubbard 
bands are present in the LDA+DMFT electronic struc- 
ture. The Fermi-liquid behavior is clear from the self- 
energy on the real-frequency axis, which we plot in Fig. [2] 
for the rfpp-hamiltonian. Although it shows a quite rich 
structure as a function of energy, the real part displays 
clear linear behavior at low-frequency. The imaginary 
part is small around u = Q and has a quadratic frequency 
dependence at low frequency. It does increase to rather 
large values at higher frequencies, however, especially for 
occupied states. Hence, our results are in general agree- 
ment with the previous calculations of Anisimov et alA 
and with the experimental photoemission (PES)^ and 
X-ray absorption (XAS)^ spectra of LaFeAsO, which 
report a moderately-correlated system with mass rcnor- 
malisation around 1.8 — 2.0. 

In order to analyze the strength of correlations for dif- 
ferent Fe 3c? orbitals we calculated the corresponding 

quasiparticle residues Z,„ = 1 — ];ni[ '^^'^^^"'' 

from the self-energy Eq. ([5]) on the Matsubara grid 
(hence, avoiding all uncertainties related to the analyti- 
cal continuation). The values are 0.609, 0.663, 0.609, and 



dncy, and degenerate dxz/dyz or- 



0.596 for the dz^ , d^^ 
bitals, resp. In this dpp energy window the Wannier func- 
tions become quite localized and their spread is expected 
to be isotropic. Indeed, within the dpp-hamiltonian the 
difference in Zm between the orbitals is rather small. 
The resulting value for the average mass renormaliza- 
tion (between 1.5 and 1.7) is in reasonable agreement 
with the experimental estimate of 1.8 extracted in Ref.[46l 
from experimental PES. The smaller mass renormaliza- 
tion found in our calculation compared to the experi- 
mental value can be attributed to the single-site approx- 
imation of DMFT. Spatial spin fluctuations, which are 
completely neglected in this approach, can eventually in- 
crease the effective mass of the quasiparticles. 

The partial densities of states for all Fe 3d orbitals 
computed within the dpp-modcl are displayed in Fig. [31 
The partial LDA-f DMFT DOS for the - and yz, xz 



FIG. 4: Momentum-resolved spectral function of LaOFeAs, 
dpp Hamiltonian. Dark areas mark large spectral weight. 



orbitals are shifted upwards relative to the xy. Indeed, 
we found that the crystal field (CF) splitting between the 
Fe 3d orbitals is somewhat affected by correlations. The 
splitting between the lowest xy and highest z^ orbitals 



remains unchanged (« 0.3 eV), while the 



y and 



yz,xz CF levels are shifted upwards by 0.15 and 0.08 eV 
relative to their positions in LDA. In LDA-I-DMFT they 
are located at 0.25 and 0.18eV, respectively, above the 
xy orbital. 

It is also instructive to look at the momentum-resolved 
spectral function A(k, w) of the crystal. It is obtained 
from the lattice Green's function, Eq. using the real- 
frequency self-energy and tracing over the orbital degrees 
of freedom. The result for the dpp-mode\ is shown in 
Fig. |4] for an energy range including Fe-d, As-p, and 
0-p states. In agreement with the Fermi-liquid picture 
of moderately-correlated quasiparticles discussed above, 
one can sec well-defined excitations around the Fermi 
level, which get more diffuse at higher binding energies. 
The bands above the Fermi level are less affected, since 
the self-energies are quite asymmetric and smaller for 
positive frequencies, see Fig. [2] Additionally, it is easy to 
see that the As-p states, dominantly in the energy range 
[—3.5, —2] eV, hybridize stronger with the Fe-d states and 
get, thus, affected by correlations. This effect is almost 
absent for the 0-p states, since they hybridize much less 
with Fe-d. 

In Fig. [5] we show a comparison between the LDA 
band structure and the LDA+DMFT k-resolved elec- 
tronic structure in a low-energy range around the Fermi 
level. This again reveals the coherent quasiparticles at 
the Fermi level, as well as more diffuse bands at higher 
energies. The crossover between long-lived quasiparticles 
and more diffuse states with a shorter lifetime is around 
— 0.4eV, in qualitative agreement with existing ARPES 
data4^. A point to mention here is the effect of the CF 
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FIG. 5; Comparison of the momentum-resolved spectral func- 
tion of LaOFeAs at low energies. Upper panel: LDA. Lower 
panel: LDA-f DMFT dpp hamiltonian. 



splitting on the band structure. For example, a differ- 
ence between the LDA and DMFT results can be seen 
for the excitation with predominantly xy character. In 
LDA it forms a holc-pockct with an excitation energy of 
-1-0.08 cV at the F point. Due to correlations, however, 
this band is shifted down significantly to the Fermi level, 
and the third hole pocket stemming from the d^y orbital 
could eventually vanish upon electron doping. 

One has to keep in mind that a direct comparison to 
experimental data is difficult for this compound, since (i) 
the experiments where done at low temperatures in the 
SDW phase, whereas our calculations arc done at room- 
temperature using the tetragonal crystal structure, and 
(ii) ARPES experiments on the 1111 family of pnictide 
superconductors are difficult to perform because of dif- 
ficulties with single-crystal synthesis. Nevertheless, on 
a qualitative level, there is a satisfactory agreement be- 
tween LDA+DMFT and experiments. 



TABLE II: Quasiparticle weights for different interaction pa- 
rameters, with the Wannier orbitals constructed for a very 
large window W = [—5.5, 13.6] eV 



Interactions 


'1 '1 '2 

z X — y xy yz,zx 


U = 3.00, J = 0.80 
f/ = 3.00, J = 0.60 
U = 3.70, J = 0.80 


0.62 0.66 0.58 0.58 
0.74 0.77 0.72 0.72 
0.58 0.61 0.52 0.56 



We also studied the dependence of the results on the 
values of the interaction parameters U and J. The re- 
sulting quasi-particle renormalizations Z,„ are listed in 
Table IT] Comparing the first two rows, one can see that 
a smaller value of J decreases the degree of correlations. 
The third line correspond to values similar to the ones 
used in Ref. 0, giving very similar results. We also in- 
creased U to the (unphysically) large value oiU ^ 5.0 eV, 
and the system still displays metallic behavior, although 
more correlated. Hence, our calculations strongly sug- 
gest that LaFeAsO is not close to a Mott metal-insulator 
transition. 

In order to check the robustness of our results, we also 
investigated the effect of increasing even further the spa- 
tial localization of the Wannier functions, correspond- 
ing to a very large energy window W = [—5.5, 13.6] eV. 
We did several calculations for different parameter sets, 
and the resulting quasi-particle renormalizations of 
all these calculations are listed in Table [III For this case, 
no cRPA calculations for the interaction matrices were 
performed, but it is expected that U and J will slightly 
increase with more localised Wannier orbitals. In that 
sense, the first row of Table HIl corresponds to interaction 
parameters that could be realised for these Wannier func- 
tions. It is very satisfying to see the calculations gave al- 
most identical quasi-particle renormalizations. Also the 
dependence on U and J is very similar to the one we 
found for the dpp hamiltonian. In that sense we consider 
our calculations to be converged in terms of the number 
of Bloch bands that are included for the construction of 
the Wannier functions and the local Hamiltonian. 



C. Remarks on calculations using the 
d-Hamiltonian and extended Wannier functions 



TABLE I: Quasiparticle weights for different interaction pa- 
rameters, with Wannier orbitals constructed from W = 
[—5.5,2.5] eV {dpp hamiltonian). The values in boldface cor- 
respond to the interaction parameters obtained from cRPA. 



Interactions 


'A '2 'I 

z x — y xy yZyZX 


U = 2.69, J = 0.79 
U = 2.69, J = 0.60 
U = 3.70, J = 0.80 
U = 5.00, J = 0.80 


0.61 0.66 0.61 0.60 

0.72 0.76 0.73 0.71 
0.52 0.57 0.53 0.52 
0.41 0.45 0.43 0.42 



In this section, we address the LDA+DMFT calcula- 
tions performed with the so-called rf-hamiltonian, where 
only the 10 Fe-d bands around the Fermi level arc used 
for the construction of the Wannier orbitals. In doing 
so, we shall shed light on the discussion which has ap- 
peared in the literaturo^i^i^i^ regarding the results of 
LDA+DMFT calculations by different authors, and the 
degree of correlations of the 1111 family of pnictide su- 
perconductors. 
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FIG. 6: (Color online) Impurity spectral function within the 
d hamiltonian using interaction parameters \J = 4.0 eV and 
J = 0.7 eV as have been used in Ref. 0. For those values a 
very correlated metal is obtained. 



clearly that the spherical approximation is highly ques- 
tionable when using only the d-bands for the Wannier 
construction. Of course, the full anisotropic interaction 
matrices can in principle be used in the LDA+DMFT 
calculation, but this raises the very delicate issue of a 
reliable orbital- dependent double-counting correction. 

Another consequence of using dclocalizcd Wannier 
functions is that they lead to significant non-local in- 
teractions Vdd, which we found to be (from cRPA) of 
order 0.23C/ to 0.32C/. These interactions are completely 
neglected in the single-site local DMFT approach, sug- 
gesting the need for a cluster extension in that case. For 
these various reasons, we have reservations against using 
a d-only Hamiltonian with extended Wannier functions 
for DMFT calculations on LaFeAsO, as also previously 
emphasized Ref. fitj 



2. Consistency with previous calculations 



1. Wannier functions and interaction matrices 

The first thing to note is that the Wannier func- 
tions constructed from a small energy window encom- 
passing only the Fe-d bands are quite extended and very 
anisotropic, as discussed in details in Ref. \^ This is di- 
rectly reflected in the interaction matrices calculated by 
cRPA in this restricted energy window: 



u: 



aa I 

mm' IcRPA 



aa I 

mm' IcRPA 



/ 0.00 
' 1.41 
1.26 
1.87 

V 1.87 

/ 3.17 
' 2.02 
1.72 
2.22 

V 2.22 



1.41 
0.00 
1.91 
1.54 
1.54 

2.02 
3.36 
2.16 
2.04 
2.04 



1.26 
1.91 
0.00 
1.33 
1.33 

1.72 
2.16 
2.17 
1.73 
1.73 



1.87 
1.54 
1.33 
0.00 
1.44 

2.22 
2.04 
1.73 
2.73 
1.84 



1.87 \ 
1.54 ' 
1.33 
1.44 
0.00 / 

2.22 \ 
2.04 ' 
1.73 
1.84 
2.73 / 



which display a strong orbital dependence. For in- 
stance, the intraorbital (Hubbard) interaction spans from 
2.17 eV to 3.36 eV. The interaction matrices in the spher- 
ical symmetric approximation using, the averages U = 
2.14 and J = 0.59, are 



( 0.00 
1.25 
1.25 
1.85 

V 1.85 

/ 2.82 
1.77 
1.77 
2.18 

V 2.18 



1.25 
0.00 
2.06 
1.45 
1.45 

1.77 
2.82 
2.31 
1.91 
1.91 



1.25 
2.06 
0.00 
1.45 
1.45 

1.77 
2.31 
2.82 
1.91 
1.91 



1.85 
1.45 
1.45 
0.00 
1.45 

2.18 
1.91 
1.91 
2.82 
1.91 



1.85 \ 
1.45 ' 
1.45 
1.45 

0.00 y 

2.18 \ 
1.91 ' 
1.91 
1.91 
2.82 / 



The largest deviation in this case is AJ7 = 0.65 eV, cor- 
responding to a relative error of about 0.26. This shows 



Nevertheless, in order to clarify apparent dis- 
crepancies between previously published LDA-I-DMFT 
results^ii^i^i^i^, we performed calculations within the d- 
hamiltonian, for several interaction parameters reported 
in the literature. For the values U = 4.0 eV and J = 
0.7 eV used in Ref. 0, we do confirm that the results 
then display very strong correlations with quasiparticle 
renormalizations ranging from Z = 0.11 (xy orbital) to 
Z — 0.34 (x^ — orbital). One may note that within 
the d-model there is a substantial orbital dependence of 
Zm, with a stronger renormalization predicted for the 
xy, yz and zx orbitals. This is a clear consequence of 
the Wannier functions being much more delocalized and 
anisotropic. 

The scattering rate at this inverse temperature of /3 = 
40eV-i is quite sizable (ImE(w+ = 0) « -0.4 ... - 0.6, 
depending on the orbital), showing that the system is 
on the verge of a coherence-incoherence crossover and a 
bad metal. The impurity spectral function is plotted in 
Fig. [6l It resembles very much the one shown in Fig. 3 
of Ref. [2, showing clear signatures of lower and upper 
Hubbard bands. There are, though, some discrepancies 
with the total weight and the positions of the Hubbard 
bands, but given the differences in the calculation (un- 
derlying electronic structure method, temperature, in- 
teraction vertex which here is only density-density) . this 
agreement with Ref. is quite satisfactory. 

Furthermore, usingthe parameters U = 0.8 eV and 
J = 0.5 eV from Ref. [j, we find renormalizations in the 
range Z k, 0.7—0.8. This is somewhat smaller than 
reported in Ref|3, although not in drastic disagreement. 

Finally, we investigated the dependence on the Hund's 
rule coupling of calculations performed with the d-only 
hamiltonian. Decreasing J to the much lower value 
J = 0.2 eV but keeping U = 4cV, we find the system 
to be much less correlated [Z between 0.63 and 0.73). 
We thus confirm, for those calculations, the great sensi- 
tivity to the Hund's coupling reported in Ref.y. We note 
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however that, although reducing J does make the system 
somewhat less correlated in this case too, this sensitivity 
is much weaker when calculations are performed with the 
full dpp hamiltonian, as reported above. 

3. Origin of the sensitivity to the Hund's coupling: level 
crossings 

In order to understand the origin of the remarkable 
sensitivity of the correlation strength to the value of J ob- 
served with the d-Hamiltonian we have studied the evo- 
lution of the ground state of the Fe 3d atomic shell as 
function of J. We obtained the 3d level positions corre- 
sponding to two different choices of the energy window: 
the " small" one corresponding to the d-Hamiltonian and 
comprising 10 Fe 3d bands and the very large one com- 
prising all As 4p, O 2p and Fe 3d bands as well as all unoc- 
cupied bands up to 13 eV above Ep. The non- interacting 
level positions e^^/ arc then obtained as 

0,0" \ ^ pa, (J a TyOL.a^ yi(T,dc 

^mni' / J ^my ^'k.v^um' mm' ' \'^'^ ) 

where the double counting term S^J^, is calculated in 
accordance with Eq. ([32|l but with the "atomic" occu- 
pancy = 6 of the Fe 3d shell. We used the same 
values of [/ = 2.14 and 2.69eV for the "small" and "very 
large" window choices, respectively, while the value of 
J was varied from 0.1 to 0.5 eV. With e^'^^, correspond- 
ing to the d-Hamiltonian we observed a level crossing at 
J w 0.2 cV with the atomic ground state changing from 
the one with spin moment S* = 1 to the one with S = 2. 
In the case of the "very large" window the ground state 
always corresponds to = 2, and the splitting between 
the ground state and first excited level is constant. It is 
obvious that a drastically different behavior of those two 
" atomic" models is related to the corresponding level po- 
sitions Cmm', which are computed using different choices 
for the Wannicr orbitals. The observed change of the Fe 
3d atomic ground state, induced by increasing J, hints 
on a possible strong dependence of correlation strength 
on the Hund's rule coupling for LaFeAsO, which is in- 
deed observed in our LDA+DMFT calculations with the 
d-Hamiltonian. However, this sensitivity stems from a 
particular choice of delocalizcd and anisotropic Wannicr 
functions and is much less pronounced when the en- 
ergy window for the Wannier function construction is 
increased. 

The bottom-line of this investigation is that all pre- 
viously published calculations seem to be technically 
correct. However, as discussed above, one introduces 
several severe approximations when dealing with the 
d-hamiltonian only, and the justification of these ap- 
proximations (restriction to local interactions, single-site 
DMFT, etc..) is questionable. This is especially true in 
this compound, due to the strong covalency between iron 
and arsenic states. 



IV. CONCLUSION AND PROSPECTS 

In the first part of this work, we present an imple- 
mentation of LDA+DMFT in the framework of the full- 
potential linearized augmented plane waves method. We 
formulate the DMFT local impurity problem in the basis 
of Wannier orbitals, while the full lattice Green's func- 
tion is written in the basis of Bloch cigenstates of the 
Kohn-Sham problem. In order to construct the Wannier 
orbitals for a given correlated shell we choose a set of lo- 
cal orbitals, which are then expanded onto the KS cigen- 
states lying within a certain energy window. In practice, 
we employ the radial solutions of the Schrodinger equa- 
tion for a given shell evaluated at the corresponding lin- 
earization energy as local orbitals. By orthonormalizing 
the obtained set of basis functions we construct a set of 
true Wannier orbitals as well as projector operator matri- 
ces relating the Bloch and Wannier basis sets. Wc derive 
explicit formulas for the projected operator matrices in 
a general FLAPW framework, which may include differ- 
ent types of augmented plain waves, lo and LO orbitals. 
Our new implementation is benchmarked using the test 
case of SrVOa , for which wc have obtained spectral and 
electronic properties in very good agreement with results 
of previous LDA-I-DMFT calculations. 

In the second part of this paper we apply this 
LDA-I-DMFT technique to LaFeAsO in order to assess 
the degree of electronic correlations in this compound 
and clarify the ongoing controversy about this issue in 
the literature. We solved the DMFT quantum impurity 
problem using a continuous-time quantum Monte Carlo 
approach. The Wannicr functions arc constructed using 
an energy window comprising Fe 3d, As 4p and O 2p. The 
resulting Wannier orbitals are rather well localized and 
isotropic. We take the average values of U =2.69 eV and 
J =0.79 from constrained RPA calculations, where the 
Wannier functions and screening channels are consistent 
with our setting of the LDA-I-DMFT scheme. We have 
checked the robustness of these results by increasing the 
size of the energy window, what resulted in a very similar 
physical picture. 

Our LDA+DMFT results indicate that LaFeAsO is a 
moderately correlated metal with an average value for 
the mass renormalization of the Fe 3d bands about 1.6. 
This value is in reasonable agreement with estimates from 
photoemission experiments. 

We also consider a smaller energy window that includes 
Fe-d states only. The resulting Wannier functions in this 
case are quite extended, leading to anisotropic and non- 
local Coulomb interactions. We take different values for 
U and J, including the ones used in previous theoretical 
LDA+DMFT approaches. We demonstrate that differ- 
ent physical pictures ranging from a strongly correlated 
compound on the verge of the metal insulator transition 
to a moderately to weakly correlated one can emerge de- 
pending, in particular, on the choice of the Hund's rule 
coupling J as observed in Ref. [^. However, there are 
conceptual difficulties when constructing a local Hamil- 
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tonian from rather delocalized Wannier orbitals. The 
interactions are very anisotropic and orbital dependent, 
and non-local interactions could also become important. 

In summary, we demonstrate that the discrepancies in 
the results of several recent theoretical works employing 
the LDA+DMFT approach stem from two main causes: 
i) the choice of parameters of the local Coulomb inter- 
action on the Fe 3d shell and ii) the degree of localiza- 
tion of the Wannier orbitals chosen to represent the Fe 3d 
states, to which many-body terms are applied. Regarding 
the first point, the calculated interaction parameters em- 
ployed in the present work arc significantly smaller than 
the values hypothesized in Refs. t^lSi]- Regarding the sec- 
ond point, we provide strong evidence that the DMFT 
approximation is more accurate and more straightfor- 
ward to implement when well-localized orbitals are con- 
structed from a large energy window encompassing Fe-3d, 
As-ip and 0-2p. This issue has fundamental implications 
for many-body calculations, such as DMFT, in a realistic 
setting. 



APPENDIX A: BENCHMARK: SRVO3 

For benchmarking purposes, we present in this ap- 
pendix LDA4-DMFT results for an oxide that has become 
a classical test compound for correlated electronic struc- 
ture calculations, namely the cubic perovskite SrVOa. 
As a paramagnetic correlated metal with intermediate 
electron-electron interactions, it is in a regime that is 
neither well described by pure LDA calculations nor 
by approaches such as LDA-I-U that are geared at or- 
dered insulating materials. From the experimental side, 
SrVOs has been characterized by different techniques 
(angle-resolved and angle- integrated photoemission spec- 
troscopy, optics, transport, thermodynamical measure- 
ments etc ) 

LDA-I-DMFT calculations have been performed both 
for an effective low energy model that comprises the three 
degenerate bands of mainly t2g character that are located 
around the Fermi energy - taking advantage of the cubic 
crystal field that singles out this group of bands - and 
for a bigger energy window comprising also the oxygen 
p-states.i^^'ia^'23^'26S^. 

In the low-energy effective t2g model a quasi-particle 
renormalization of Z ~ 0.6, compatible with experi- 
ments, is obtained for U values around 4eV. The re- 
maining spectral weight is shifted toward lower and up- 
per Hubbard bands. The lower Hubbard band, located 
around — 1.5eV binding energy has indeed been observed 
in photoemission; the high energy satellite of the 
model is located around 2.5 eVi^ 

Concerning calculations taking into account also the 
ligand states, it should be noted that possible LDA errors 
on the separation of p- and d-states are not corrected by 
DMFT, since only the d-states are treated as correlated. 

In the present work, we use SrVOa as a benchmark for 
our projector orbitals implementation of LDA-I-DMFT, 




CO /eV 

FIG. 7: (Colore online) DOS for SrVOa, d-only model (small 
energy window). Top panel (a): Total DOS of LDA (black) 
and LDA-I-DMFT (red). Bottom panel (b): LDA local or- 
bitals (black), impurity spectral function A{u}) (red), and 
vanadium t2g partial DOS (green). Coulomb parameters for 
these calculations are given as inset. 
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FIG. 8: (Color online) Total DOS for SrVOs, dp Hamiltonian 
(vanadium t2g and oxygen p). 

with results very similar to previous theoretical studies. 
We performed two kinds of calculations: (i) We used as 
an energy window the range from — 1.35eV to 1.90 eV 
which comprises the t2g bands, and - at some k-points - 
one or both of the Cg bands. This is closest in spirit to a 
t2g model within a Wannier function formalism, though 
not exactly the same due to the inclusion of some Cg con- 
tribution. To recover a Wannier prescription one would 
in fact have to choose a k-dependent window, such as to 
include exactly three bands at each k-point, correspond- 
ing to the three-fold degenerate manifold of dominantly 
t2g bands, (ii) We used an energy window of — S.lOeV to 
1.90 eV, spanning both, the bands used in (i) and the oxy- 
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gen p dominated bands located between —8 and — 2eV. 

Please note that, in order to be consistent with existing 
literature, we use a different parametrisation of the inter- 
action matrix compared to Sect. lIIII Here we define U to 
be the onsite intraorbital Coulomb interaction, U—2J to 
be the inter-orbital interaction for electrons with oppo- 
site spin, and U—3J the inter-orbital interaction between 
electrons with equal spin. 

The results for the first case are shown in Fig. [T] 
The upper panel. Fig. [7] (a), displays the total spec- 
tral function of the thus defined model within LDA and 
LDA-I-DMFT. Our results recover previously published 
results, with a quasiparticle renormalisation of around 
Z = 0.60 for a values U = 4.0 eV and J = 0.65 eV. The 
contribution of the Cg bands to the total DOS can easily 
be identified from the LDA-f-DMFT spectra, where an 
additional hump between the quasiparticle peak and the 
upper Hubbard band appears. 

Fig. [7] (b) shows the local orbitals used for the DMFT 
calculations, together with the corresponding impurity 
spectral function A{ll}) and the vanadium t2g partial 
DOS. The latter one is obtained bu projecting the lattice 
Green's function to character using the partial pro- 
jectors to be introduced in App. |B] The main difference 
to panel (a) is the absence of the additional Cg character. 

Finally, Fig. M shows the LDA-hDMFT spectral func- 
tion compared to the LDA density of states, as calculated 
within the larger energy defined in (ii) above. Since the 
Wannier function are more localised, as compared to case 
(i), the value for the Coulomb interactions has to be ad- 
justed accordingly, and we chose a value oi U = 6.0 eV. 
As can be seen in the figure, ligand states are barely 
modified by the correlations, and the results for the t2g- 
derived bands are very close to what is seen in the effec- 
tive low energy model. The quasiparticle renormalisation 
is Z = 0.57, in good agreement with the pure t2g treat- 
ment discussed before. 

The results of these calculations correspond to what 
can be expected on the basis of previously published 
work, and thus validate our new implementation. 



inside the sphere uf{r, Eii)Yj^x<y, Eii)Y^x<y and 

uf {r, Ei2)Yj;^Xcr which do not form an orthonormalized 
basis set. It is more convenient to construct these pro- 
jectors if the wave function is rewritten in an orthonormal 
basis set. 

In a general form, inside a given sphere we can express 
V^L(r) as 

V'L(r) = + ^i'™"' + ^'™"'2, (Bl) 

Im Im Im 

where we simplify the notation by omitting the angu- 
lar and spin parts and defining A'l^, and C;',„ as 
combined coefficients which are generally k-dependent 
and contain the sum over the plane waves and local or- 
bitals. We also define un =ui[r, En), iii = ui{r, En) and 
ui2 = ui{r,Ei2). 

We then rewrite ipi^^i'r) as a function of a set of or- 
thogonal orbitals 4>j{r), j=l,2,3 as follows 

Im j 

(B2) 

The coefficients d™ are the matrix elements of the 
square root of the corresponding overlap matrix 

/ 1 {un\ui2)\^ 

C= {iiil ill) {ui\ U12) . (B3) 

\ {Ul2\ Uii) (U/2I Ul) (U/2I U12) I 

In this way, rewriting Eq. IB2I as 

ruv) = YY~^TM'-)- (B4) 

Im j 

the matrix elements of the projector to a given atom with 
Im character finally reads. 



APPENDIX B: 



PROJECTORS FOR PARTIAL 
DOS 



In order to calculate the partial density of states for a 
given atomic site and particular orbital character (corre- 
lated or not) we construct a different type of projectors, 
which we call 0*''^. 

The Wannier operators of Eq. [M] project onto a given 
Wannier-like orbital. On the other hand, the new set 
0*''^, as we will show, project onto a given orbital of cer- 
tain character for which we do not apply any orthonor- 
malization process as in the first. Unlike the Wannier 
projectors, the ©''"^'s can also project to other orbitals 
atoms apart from the correlated set. 

A given orbital character contributes in the eigen- 
states through the solutions of the Schrodinger equation 



~lm 



(B5) 



The spectral function of a given atom i with orbital 
character m, is obtained as 



A^^{]^,u;)^—Im 

TT 



E©-'.WG-,,(k,a;+)e:;;^,^(k) 



APPENDIX C: INFLUENCE OF THE 
ROTATIONAL INVARIANCE OF HUNDS-RULE 
COUPLING IN MULTI-ORBITAL SYSTEMS 

In our DMFT calculations using CTQMC as impu- 
rity solver, we restricted the Hund's rule interaction to 
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FIG. 9: Quasi-particle renormalisation Z in a 3-orbital Hub- 
bard model. Calculations have been done using a semicircular 
DOS with bandwidth W = 4i. Top panel: n = 3 (half filling). 
Bottom panel: n — 2. 



Ising-type interactions only, although there is no concep- 
tual limitation of the algorithm to this type of interac- 
tions. The reason for doing this is of purely technical 
nature, since in this case one can diagonalize the local 
problem very efhcicntly, and furthermore, it enables us to 
use the so-called segment-picture update scheme,— which 
increases the efficiency of the CTQMC method a lot. 

One may now ask how results change if the fully 
rotational-invariant Hund's rule exchange is taken into 
account. For this purpose, we study a multiband model 
Hamiltonian, assuming degenerate bands, no interband 
hybridizations, and a semicircular density of states. Ap- 
plying the self-energy functional theory (SFT)f^ we can 
study the quasiparticle renormalization Z as function of 
interactions U and J. In this study, we choose the con- 
vention of setting the intraorbital Coulomb repulsion to 
U and the interorbital to U' = U — 2 J, and give all en- 
ergies in units of the single-particle hopping amplitude t, 
i.e., the band with of the DOS is = U. 

In addition to the density-density interactions, we con- 
sider also the additional spin-flip and pair-hopping terms 
of the local Hamiltonian, 



J 



(CI) 



(C2) 



We do calculations at T = 0, and choose the reference 
system for the SET framework to consist of one bath de- 
gree of freedom for each correlated orbital. Hence, going 
up to M = 5 orbitals, we have to diagonalize a local 
problem consisting of at most 10 orbitals. 

The upper panel of Fig.[5]shows Z for a 3-orbital model 
at half-filling, n = 3, for J = O.IU. A tremendous re- 
duction of the critical Uc of the metal-to-insulator tran- 
sition (MIT) is observed, already for Ising-like interac- 
tions. This is a well known fact that for multi-orbital 
systems at or close to half-filling, the effect of J should 
be strongest i^Si^iS The inclusion of spin-fiip and pair- 
hopping terms gives raise to two effects, (i) For moderate 
correlations, Z « 0.6, these terms lead to a slight reduc- 
tion of Z, but (ii) the critical U for the MIT is shifted 
upwards. This qualitatively holds also away from half- 
filling, which can be seen in the lower panel of Fig. [HI 
where we plotted Z for ?? = 2. Although the transition 
is not of first order any more, one can again identify 
two regimes. For moderate correlations, Z decreases, 
whereas close to the transition the spin-flip and pair- 
hopping terms increase the renormalization Z and the 
critical Uc is pushed to higher values. This is consistent 
with a numerical renormalization group study for the 2- 
orbital Hubbard modeli^ 

We also considered the case relevant for pnicitde ma- 
terials, i.e., M = 5, n = 6, and Z around 0.5. This 
regime could be realize by setting (i) U ~ 3.5t and 
J = 0.35^, which shows a reduction of Z from 0.52 to 
0.47 due to spin-flip and pair-hopping, or (ii) U = 2t, 
J = OAt, giving a reduction from 0.61 to 0.57. In conclu- 
sion, this analysis shows that the picture of a moderately- 
correlated metal as argued in Sect. IIIII holds also when a 
fully rotational-invariant Hund's exchange is considered. 
For systems close to a MIT, this is no longer true and 
the spin-flip and pair-hopping terms become crucial. 
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